Subject classification: 77. 84. Fa 
Institute for Condensed Matter Physics 
National Academy of Sciences of Ukraine 
Uzhgorod State University 
Physics Departmental 

Hydrostatic pressure influence on dielectric permittivity of KH2PO4 and KD2PO4 
in the piezoelectric resonance region 

By 

R. R. Levitskii 1 , A. G. Slivka 2 , A. P. Moina 1 , P. M. Lukach 2 , A. M. Guivan 2 

We measured transverse and longitudinal dielectric permittivities of KH2PO4 and KD2PO4 in the piezo- 
electric resonance frequency region under the hydrostatic pressure. The transverse permittivity decreases 
with pressure in the paraelectric phase and increases in the ferroelectric phase. The pressure dependence 
of the transverse permittivity of deuterated KD2PO4 is well described by the presented theory. The longi- 
tudinal permittivity of a pure KH2PO4 exhibitsseveral peaks in the vicinity of the transition point. Upon 
increasing hydrostatic pressure, the peaks get closer and may merge. 



1 Introduction 

An important role in the phase transition and dielectric response of ferroelectrics and antiferro- 
electrics of the KH2PO4 family is played by the geometry of hydrogen bonds in these crystals. 
Pressure dependence of the separation S between equilibrium proton (deuteron) sites on a bond 
governs the corresponding dependences of the transition temperature and longitudinal static di- 
electric characteristics of the crystals [Q, @, ||. Moreover, the dependence Tq on S is universal || || 
for several crystals of the KH2PO4 type having a three dimensional network of hydrogen bonds 
(MeD 2 X0 4 , Me = K, Rb, ND 4 , X = P, As, KH 2 P0 4 , and NH 4 H 2 P0 4 ). 

High pressure studies are one of the most useful methods of exploring the role of hydrogen 
subsystem geometry in the phase transition and dielectric response of these crystals. Influence of 
hydrostatic pressure on the transition temperature, spontaneous polarization, and Curie constant 
of the KH2PO4 family ferroelectrics is well studied experimentally. This influence for the prototype 
compounds KH2PO4 and KD2PO4 was theoretically described yet in Rcfs. |l], || using the four- 
particle cluster approximation for the proton ordering model. The fact that this model can account 
for the external pressure effects is considered as one of its experimental evidences. 

However, much less is known about pressure effects on the transverse dielectric characteristics 
of these crystals. Only in Ref.j^] the variation of transverse permittivity e a of an undeuterated 
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KH2PO4 with hydrostatic pressure at room temperature was measured. It was established that e a 
decreases with pressure, and the slope ds a /dp decreases as frequency is increased from 400 Hz up 
to 25 kHz. As far as we know, there is no experimental data for influence of external pressures on 
transverse dielectric properties of deuterated KH2PO4 type crystals. Neither it is known in what 
way external pressure affects the longitudinal or transverse dynamical dielectric characteristics of 
these crystals. 

The aim of the present study is to fill this gap, namely, to explore the dependences of transverse 
and longitudinal dynamic dielectric permittivities of the KH2PO4 and KD2PO4 on hydrostatic 
pressure in a wide temperature range. Measurements are carried out at the frequency of IMHz, 
that belongs to the piezoelectric resonance frequency region of these crystals H . This fact should 
not affect the behavior of the transverse permittivity, provided the samples are oriented precisely, 
but the temperature curves of the longitudinal permittivity must be essentially different from the 
static (0 - 10 4 Hz) and high-frequency (I0 9 - 10 12 Hz) ones. The pressure dependence of the 
transverse permittivity of a deuterated KD2PO4 is described within a microscopic theory. Thence, 
we can find out whether the pressure changes in the hydrogen bonds geometry determine the 
pressure dependence of not only the longitudinal |§[ gj, but also of the transverse dielectric 
response of this crystal. 



2 Experimental methodics 

We measured the dielectric permittivity of two crystals: a pure KH2PO4 (transition temperature 
at ambient pressure Tqo ~ 123 K) and a highly deuterated sample with Tea ~ 210 K (0.87 
nominal deuteration) , hereafter abbreviated as KD2PO4. The obtained decrease of the transition 
temperatures with hydrostatic pressure dTc/dp is —4.6 K/kbar in KH2PO4 and —2.1 K/kbar in 
KD2PO4, that well accords with the literature data J7j. 

The permittivity e was determined from the samples capacity using the well-known formula 

£ ~e S C > 

where d is the sample thickness; C is its electric capacity; S is the area of electric contacts; 
£0 = 8.85 • I0~ 12 F/m. Capacity was measured by the conventional bridge method using the ac 
bridges E7-I2 (working frequency 1 MHz) and P5016 (10 kHz and 50kHz) within 0.2-0.4%. As 
electric contacts a silver paste was used. 

Optic grade samples were placed in a autonomous hydrostatic pressure chamber, with silicone oil 
serving as a pressure transmitting medium. Pressure was measured by mechanical and manganin 
manometers with an accuracy ±2 MPa. Temperature was measured to ±0.1 K with a copper- 
constantan thermocouple. Capacities were obtained in a dynamical regime with a temperature 
change rate 2 • I0 _2 K/c. 



3 Experimental results 

Closeness of the measuring frequency to the piezoelectric resonance region, as expected, does not 
essentially influences the transverse permittivities of the KH2PO4 and KD2PO4 ferroelectrics. 
The measured curves have a typical form of the static transverse permittivity in these crystals (see 
fig. [0). In the paraelectric phase £n slowly increases on lowering temperature, reaching its maximal 
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value at T « Tc + 15 K, somewhat decreases in the interval Tq < T < Tq + 15 K, has a jump at 
the transition and gradually decreases to a certain limiting value as temperature is lowered down 
in the ferroelectric phase. This gradual decrease in a deuterated KD2PO4 is much faster than in 
a pure KH2PO4, reflecting the fact that in KH2PO4 the phase transition is close to the second 
order, while in KD2PO4 a pronounced first order phase transition takes place with a significant 
jump of spontaneous polarization M. 

The obtained pressure dependences of the transverse permittivities of ferroelectric KH2PO4 
and KD2PO4 are analogous to the corresponding dependences of the transverse permittivity in 
antiferroelectric NH4H2ASO4 and ND4D2ASO4 M. Under hydrostatic pressure, the magnitudes of 
transverse permittivities of both crystals at temperatures far above the transition points decrease 
with dlnsn/dp = -S.G^okbar- 1 in KH 2 P0 4 at T = 160 K and -LWokbar" 1 in KD 2 P0 4 at 
T = 260 K. In the ferroelectric phase en increases with pressure, but this effect results from the 
changes in the distance of a given temperature point from the transition temperature, whereas the 
low temperature limit of en does not depend on pressure. 

In figs. and ^| we presented the temperature curves of the longitudinal dielectric permittivity 
of a pure KH2PO4 and undeuterated KD2PO4, respectively, at different values of hydrostatic 
pressure. Note an interesting behavior of the permittivity of KH2PO4, which exhibits three clear 
maxima (correspondingly, two minima) in the vicinity of the transition point. Possibly due to the 
narrowness of the temperature interval in which the peaks of the longitudinal permittivity occur, in 
a deuterated KD2PO4 the three-maximum structure was detected only at p = 1.9 kbar (see fig. ||). 
For the other pressures, only one minimum of the permittivity was clearly observed. However, 
the bends in the curves of e33 in the temperature interval between the detected paraelectric and 
ferroelectric peaks indicate that the second, ferroelectric phase minimum must occur in KD2PO4 
as well. 

As follows from fig. [|a, where the curves of longitudinal and transverse permittivities e33 and 
en are compared, the right maximum occurs at the transition point. The two other maxima and 
the left minimum, which take place in the ferroelectric phase, are reproduced at different values of 
hydrostatic pressure as well as at heating or cooling (see fig. [|b). At lower frequencies (v = 10 kHz 
and v = 50 kHz), the e33 of deuterated KD2PO4 has a typical form of the static permittivity with 
one peak at the transition point and the Curie- Weiss behavior in the paraelectric phase (fig. ^3). 

Under hydrostatic pressure, the magnitude of the transition point maximum in KH2PO4 does 
not change, while that of the left ferroelectric one is reduced from ~ 770 at lObar down to ~ 400 at 
p = 2 kbar. Beside, the temperature interval where the extrema of the permittivity occur becomes 
narrower, so that the two-minimum structure of e33 transforms to a single-minimum one. 

At temperatures below the peaks, the longitudinal permittivity of deuterated KD2PO4 has a 
typical shoulder-like form, caused by the multi-domain structure of the samples. Qualitatively 
similar dependences can be traced in a pure KH2PO4 as well, but the "fine structure" of e33 
here is fancier. Hydrostatic pressure almost does not affect the shape of the temperature curve 
£33(T < Tc). Width of the plateau is about 30 K in both crystals, with the height of about 400 
being practically pressure independent. 

As seen in figs. |^ and|3| the Curie- Weiss law is obeyed for e 33 1 of KH2PO4 and KD2PO4 above 
the transition point in a rather wide temperature range, except for a vicinity of the transition point. 
The Curie constants are very close to the corresponding static ones and decrease with pressure 
linearly from C — 2949 K in KH2PO4 and C = 3770 K in KH2PO4 at ambient pressure with the 
slopes d\nC/dp = -OTS^okbar" 1 in KH 2 P0 4 and -1.5%kbar -1 in KD 2 P0 4 (the literature data 
for the Curie constants of static permittivities are C = 2910 K and d\nC/dp = — 0.66%kbar _1 in 
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KH 2 P0 4 and 3700 K and -l^^okbar- 1 in KD 2 P0 4 , respectively §). 

Note a strong dependence of the logarithmic pressure derivatives of the dielectric characteristics 
of the crystals on deuteration and the fact that this dependence is different for the transverse 
and longitudinal quantities: \d\neu/dp\ in KH2PO4 is much higher than in KD2PO4, whereas 
for \dh\C / dp\ the reverse holds. A similar effect has been revealed earlier |t]] for the transition 
temperature, saturation polarization and Curie constant of these crystals. 



4 Theory 

We restrict our theoretical calculations by the case of highly deuterated KD2PO4. The known mod- 
els of dielectric relaxation in hydrogen-bonded crystals JlO], jn|, based on Glauber's Gj| dynamics 
of pseudospins, are suitable for description of the high-frequency permittivities (10 — 10 12 Hz) 
only, when a crystal is effectively clamped [jl2]. Therefore, in order to describe the above presented 
pressure dependences of the longitudinal dielectric permittivity one must develop a special model 
of dielectric relaxation which would take into account dynamics of piezoelectric strain Sq as well. 
However, since this dynamics is not essential for the transverse dielectric response of KD2PO4, the 
latter can be described within static models. We shall do that within the proton ordering model 
modified to the case of strained crystals || . 

Calculations are performed with the conventional four-particle cluster Hamiltonian 
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Two eigenvalues of the Ising pseudospin a q f = ±1 are assigned to two equilibrium deuteron sites 
on the f-th bond in the q-primitive cell. The internal fields z q f include a long-range interactions 
between deuterons (dipole-dipole and lattice mediated) taken into account in the mean field ap- 
proximation, cluster fields A q f which describe the short-range interactions with the pseudospin a q f 
not explicitly included into the cluster Hamiltonian, and an external electric field applied along 
the a-axis of the tetragonal unit cell 
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fj} q j is the projection on the a-axis of the effective dipole moment created by displacements of heavy 
ions and redistribution of electron density induced by deuteron ordering. In the case of transverse 
electric field E±, the following relations between the mean values of pseudospins and between the 
effective dipole moments /xjj are obeyed 

(c 9 2) = (094), Z24 = Zq2 



V24 

Ml = Mjl = _ M ? 3) Mq2 = M^4 = 0- 
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Parameters of the short-range deuteron correlations U, V, $ are functions of the so-called Slater 
energies e, w, W\: 
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Pressure dependences of the Slater energies and parameters of the long-range interactions are 
modelled as H 
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and 



J ff ,( qq ') = J$(qq') 
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(5) 



Here we take into account the quadratic dependence of these parameters on the separation S 
between two equilibrium deuteron sites on a bond, whereas the pressure dependence of (5 is known 
to be linear [fl4| 

2 :! 



5 = S + S lP = S ( 1- ^^ej ) , 



3 = 1 



5 = 12iy=i ^i?- According to ^ [| we suppose that the only essential mechanism of the pressure 
influence on the Slater energies is a decrease in the D-site distance 6. However, for the parameters of 
the long-range interactions there exist other important mechanisms of variation with pressure (for 
instance, the dipole-dipole interactions increase when the distance between deuterons is reduced), 
taken into account in (|^) via the expansions in the diagonal components of the strain tensor 
(* = 1,2,3). 

Excluding the fields A q f from Zqt by making use of the condition that the mean values rj^ = 
{(T q f) calculated with Hamiltonian (nt) and with the one-particle Hamiltonian 
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coincide, we obtain 
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where the following notations are used 
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J%j(0) is the long-range interaction matrix Fourier transform. 

From the field-induced polarization of the crystal along the et-axis 

P = yjjf M _ ^M + rg) 



the transverse static dielectric permittivity ef(0,T,p) of a clamped crystal (at =const) can be 
calculated 
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where 
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eioo is a high-frequency contribution to the permittivity. The dielectric permittivity of a free 
crystal (at p=const) ef (0,T, p) is related to (0) by 



ef 1 (0,T,p) = e s n {0,T,p)+A7rd 2 u c^, 

where d±4 is the piezomodule, and cf^ is the elastic constant of a short-circuited crystal. 
For the involved quantities, the following relations hold in the zero external field limit 
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The order parameter r\ is found by minimization of the free energy 
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whereas to find the strains we should solve the system of equations 



3=1 



(9) 



where Cy are the elastic constants of the whole crystal, being determined from an experiment. Con- 
tribution of a deuteron subsystem to pressure or temperature dependences of the elastic constants 
is neglected. 
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5 Numerical analysis 



The values of the theory parameters for a deuterated crystal KD2PO4 with the transition temper- 
ature at ambient pressure Xco = 210 K are presented in Tables [IJ||. They were found in Ref. || 
and used for a description of the uniaxial pressure p — — 03 dependences of the transition temper- 
ature, longitudinal dielectric permittivity, and spontaneous polarization of the crystal. The clastic 
constants of the paraelectric crystals coincide with the experimental data of Ref. jl5) . The new 
parameters v a , f± — (a*i) 2 / w i an d £ioo are found by fitting the theoretical temperature curve of 
the transverse permittivity to experimental data (y a sets the slope den/dT, whereas /1, and eioo 
give the magnitude of e\\. 

Table 1: The theory parameters for a KD2PO4 crystal with T C o = 210 K, dT c /dp = -2.1 K/kbar. Plus 
and minus indices denote the quantities used in the paraelectric and ferroelectric phases, respectively. 
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The so-called deformation potentials ipai, which enter only the expression for the transverse 
permittivity, do not essentially affect the permittivity. Therefore, for the sake of simplicity we put 
ipai — 0. Contributions of the double-ionized deuteron configurations (with four deuterons close to 
a given PO4 group or with none) are neglected as well by taking w\ — > 00. 

Hence, the only new theory parameter, governing the pressure dependence of the calculated 
characteristics, is the ratio 61/60 - the relative rate of the pressure changes in the D-site separation 
6. We choose it so that at all other parameters unchanged, the correct dependence of the transition 
temperature on the hydrostatic pressure dTc/dp = —2.1 K/kbar is obtained. 

In Ref. § we have shown that at the chosen analogously values of 61/60, the recalculated depen- 
dences Tq (5) for six deuterated ferroelectric and antiferroelectric crystals with a three dimensional 
network of hydrogen bonds MeD2X04 (Me = K, Rb, ND4, X = P, As) form a single universal 
linear dependence, while Tq and 5 can be altered by either hydrostatic or uniaxial p — —03 (for 
KD2PO4) pressures. Experimentally this universality has been established by R. Nelmes et al 
[g for undeuterated KH 2 P0 4 and NH 4 H 2 P0 4 and deuterated KD 2 P0 4 and ND 4 D 2 P0 4 . At the 
adopted in this paper value of 61/60 for KD 2 PC>4 with Tco = 210 K, the dependence Tc(p)[6(p)] 
for this specific crystal also accords with the universal line [[| . 

The slopes 8^/ dp can be determined without introducing into the theory any extra fitting 
parameter on the basis of the following simple half-empirical speculations, nevertheless yielding a 
fair agreement with the experiment. 

It is believed that the deuteron ordering in the system results in displacements of heavy ions and 
electron density which contribute to crystal polarization. Since, when ordered, a deuteron shifts 

Table 2: Elastic constants of the considered crystal (units of 10 5 bar). 



c ll c 12 c 13 c 33 c ll c 12 c 13 c 22 c 23 C 33 

6.93 -0.78 1.22 5.45 6.8 -0.78 1.0 6.99 1.0 5.3 
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from its central position on a hydrogen bond to the off-central one by a distance 6/2, according 
to || we assume that the heavy ions displacements are also proportional to 6 and that [ii is 
proportional to the corresponding lattice constant a% . This yields 



1 dfij _ <5i 
M° dp 6 



(10) 



The pressure dependence of the effective longitudinal dipole moment fi^, calculated with ( |Io| ) 
provides a fair agreement with the available experimental data for the pressure dependence of 
spontaneous polarization of KD2PO4 and of the static dielectric permittivities of KD2PO4 and 
RbD 2 P0 4 j|. 

In fig. ||we depict the theoretical curves of the transverse static dielectric permittivity of 
KD2PO4 along with the experimental points of Section ||. The permittivity is calculated with (^), 
since due to smallness of the piezomodule d\4 (l6) the difference between transverse permittivities 
of clamped and free KD2PO4 can be neglected; the pressure dependence of the dipole moment /ii 



is given by (10) 



The theory qualitatively well reproduces the temperature curve of the permittivity, including 
a decrease with temperature in the paraelectric phase and the jump at the transition point, but 
is not able to explain the existence of a small broad maximum of en at temperature slightly 
higher than the transition point. That is a common drawback of all existing calculations of the 
transverse permittivity, that can be removed by assuming a temperature dependent high-frequency 
contribution to the permittivity Exoo- 

As one can see, a satisfactory quantitative agreement with experimental data for the rates of 
a decrease in the permittivity with pressure in the high-temperature phase and of an increase in 
the low-temperature phase is obtained. The fact that the value of the ratio 61 /6q, yielding the 
correct theoretical pressure dependence of transition temperature, provides the correct pressure 
dependence of the transverse dielectric permittivity, no extra fitting parameters being introduced 
into the theory, gives one more evidence on the important role played by the hydrogen bonds 
geometry, namely the separation 6 between equilibrium deuteron sites, in the phase transition and 
dielectric response of the hydrogen bonded crystals. 



6 Concluding remarks 

We performed experimental studies of hydrostatic pressure influence on transverse and longitudinal 
dielectric permittivities of ferroelectric KH2PO4 and KD2PO4 crystals at 1MHz, the frequency that 
belongs to the piezoelectric resonance region. 

The obtained temperature curves of transverse permittivity en are similar to the corresponding 
curves of static permittivity of the crystals, while the pressure dependences of the permittivities 
are analogous to those in the antiferroelectric crystals NH4H2ASO4 and ND4D2ASO4. 

Unlike en, the temperature curves of the longitudinal permittivity 833 at 1MHz are qualitatively 
different from static ones. The multipeak structure of the permittivity in the vicinity of the 
transition point is observed. Under hydrostatic pressure, the two minima of the longitudinal 
permittivity are getting closer, so that its two-minimum structure transforms to a single-minimum 
one. Above the transition point, the longitudinal permittivities of KH2PO4 and KD2PO4 obey the 
Curie- Weiss law; the pressure dependences of the corresponding Curie constants well accord with 
the data for the Curie constants of static permittivities. 



8 



The measured rates of decrease with pressure in the transition temperature and permittivity 
en in the paraelectric phase and of a increase in en in the ferroelectric phase in a deuterated 
KD2PO4 are well described within the proton ordering model. Theoretical pressure dependences 
of Tc and £11 to a great extent are governed by the ratio Si /So- This parameter denotes relative 
pressure changes in the separation S between equilibrium deuteron sites on a bond. Obtained in the 
present paper theoretical description of the pressure dependence of en, a similar agreement with 
experimental data for spontaneous polarization and longitudinal static dielectric permittivity H, 
as well as the previously revealed universality of ^, |J of the transition temperature vs H(D)-site 
distance 8 in the MeD 2 X0 4 (Me = K, Rb, ND 4 , X = P, As), KH 2 P0 4 and NH 4 H 2 P0 4 crystals, 
indicate that an important role in the phase transition and dielectric response of the crystals is 
played by the geometry of hydrogen bonds, in particular, by the separation 5 between hydrogen 
sites. 

So far we have no clear explanation of the multipeak structure of the longitudinal permittivity. 
Apparently, it is a dynamic effect connected with a piezoelectric resonance phenomena, since at 
lower frequencies as well as at frequencies above the piezoelectric resonance but below the dielectric 
dispersion region [fl7| the £33 has a typical form of the permittivity with one peak at the transition 
point and the Curie- Weiss behavior in the paraelectric phase. 

The maxima of £33 in the ferroelectric phase might be merely the peaks at different piezo- 
electric resonance frequencies (harmonics). The harmonics frequencies are determined by sample 
dimensions and by the appropriate elastic constants (only by c§ 6 for a 45° Z cut sample). The cjf 6 
has a peculiarity at the phase transition point dropping from about 6 • 10 10 dyn/cm 2 to zero 
and increasing back to this value in the paraelectric phase. Therefore, at changing temperature 
in the vicinity of the transition point, the resonant harmonics frequencies span a wide frequency 
range, so that at certain temperatures they may coincide with the measuring frequency 1 MHz. 
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Figure 1: Transverse dielectric permittivities of KH2PO4 and KD2PO4 as functions of temperature at 
different values of hydrostatic pressure. The frequency is v = 1MHz. Symbols are experimental points, 
lines are drawn for clarity. 



10 




100 110 120 130 



T(K) 



500 


p-l kbar 


\ 0.012 


400 




? 1 0.008- 


300 






200 






100 







600 
500 
400 
300 
200 
100 



p-0.5 kbar 


la 0016 




jil 0.012- 








S \° 08 " 






¥ fc 33 





90 100 110 120 130 

J 0.012- 

0.0C 




90 100 110 120 130 



T(K) 



90 100 110 120 130 T(K) 



Figure 2: Longitudinal dielectric permittivity of KH2PO4 at v = 1MHz (cooling) as a function of tem- 
perature at different values of hydrostatic pressure. Symbols are experimental points, lines are drawn for 
clarity. 
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Figure 3: Longitudinal dielectric permittivity of KD2PO4 at v — 1MHz (heating) as a function of tem- 
perature at different values of hydrostatic pressure. Symbols are experimental points, lines are drawn for 
clarity. The insets show the permittivity in the vicinity of the transition point. 
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Figure 4: Dielectric permittivities of KH2PO4 and KD2PO4 as functions of temperature at ambient 
pressure. Measuring frequencies and direction of temperature changes (heating or cooling) are indicated 
at figures. Symbols are experimental points, lines are drawn for clarity. 
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Figure 5: Transverse dielectric permittivity of KD2PO4 as a function of temperature at different values 
of hydrostatic pressure (kbar) 1 - 0; 2 - 4.1. Lines are calculated theoretically; symbols are experimental 
points of the present work. 
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